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We study the surface state of a doped topological crystalline insulator in the superconducting 
state. Motivated by Sni-a^IuccTe, we consider fully-gapped pair potentials and calculate the sur¬ 
face spectral function. It is found that mirror-protected zero-energy surface Andreev bound states 
(SABSs) appear at the (001) surface. The gapless points of these SABSs appear on the mirror 
symmetric line on the surface Brillouin zone while the positions of the gapless points depend on the 
chemical potential. In addition, due to the presence of the Dirac surface states in the normal state, 
the dispersion of the SABSs drastically changes with the chemical potential. 


I. INTRODUCTION 

Surface Andreev bound state (SABS) is an impor¬ 
tant physics in superconductors. [l[-Q. SABSs appear on 
the surface of an unconventional superconductor with a 
pair potential that changes its sign on the Fermi sur¬ 
face. 0,1113 These SABSs have been detected by tun¬ 
neling spectroscopy measurements in cuprate [il, m and 
Sr 2 Ru 04 [8|,[9|. The topological origin of SABSs has been 
clarified based on the topological invariants defined in 
the bulk Hamiltonian. Because of bulk-boundary corre¬ 
spondence, superconductors with gapless SABSs are now 
called topological superconductors [l0|, UjJ , and the con¬ 
cept of topolog ical superconductivity has been expanding 
widely [ifilf. 

On the other hand, topological insulator has been of 
great interest because of its novel surface state protected 
by time-reversal symmetry [2314^ . The superconducting 
topological insulator (STI), which is a superconductor re¬ 
alized in a doped topological insulator, has been also at¬ 
tracting a great research interest following the first obser¬ 
vation of zero-bias conductance peak that indicates the 
existence of SABSs in Cua,Bi 2 Se 3 [26|. And now, deter¬ 
mination of the pairing symmetry of Cua,Bi 2 Se 3 is recog¬ 
nized as an important problem |27l-l4^. In the STI, the¬ 
oretical studies have revealed that the cone shape SABS 
becomes twisted dispersion due to the mixing with the 
surface state in the normal state The struc¬ 

tural change of the dispersion is considered as an im¬ 
portant property for the interpretation of the tunneling 
spectroscopy. The existence of peculiar SABSs has been 
also predicted in the doped topological materials such as 
Dirac semimetal and Weyl semimetal jdTl-E^. 

Recently, a new type of topological material, topolog¬ 
ical crystalline insulator (TCI), has been proposed |^ . 
TCIs are the materials that have surface states protected 
by point group symmetry, e.g. mirror reflection sym¬ 
metry and rotational symmetries. Up to now, SnTe, 
Pbi_a,Sna,Te and Pbi_a,Sna,Se heve been revealed to be 
the TCI with non zero mirror Chern number [bll-lb^. 


They host double Dirac cones on the (001) surface pro¬ 
tected by mirror reflection symmetry [see Figdjc)]. Since 
the surface states of the TCIs are protected by point 
group symmetry, the Dirac cone can be controlled by the 
electronic field or strain as well as the magnetic field. The 
additional tunability of the surface state can be applica- 
ble to the high-speed topological transistor |55l-[^. For 
TCIs, superconductivity has been reported in In-doped 
SnTe [m . Point contact measurements in superconduct¬ 
ing Sni_a,Ina,Te show zero-bias conductance peak that 
may originate from SABSs [59|. In other words, this 
material is a candidate of a topological superconductor. 
Furthermore, it is expected that superconducting TCIs 
(STCIs) exhibit a new type of SABS since the TCIs host 
the unique surface Dirac cone in the normal state. 

Using the similarity of low energy electronic state to 
Cua,Bi 2 Se 3 , Sasaki et. al. have investigated the supercon¬ 
ducting state of Sni_a,Ina,Te with an effective model for 
around each L point (see Fig. [T]) [59|. However, dif¬ 
ferent from Cua,Bi 2 Se 3 which has a single Fermi surface, 
Sni_a,Ina,Te has four hole pockets at around each L point. 
Then, two hole pockets are projected to the same posi¬ 
tion in the surface Brillouin zone when we see the (001) 
surface. Therefore, it is inevitable to consider the over¬ 
lapping of the Fermi surfaces to reveal the surface state 
properly. 

In this paper, we study the surface states of the 
STCI using a tight-binding model which can describe 
the electronic states for the entire Brillouin zone. We 
introduce possible fully gapped pair potentials for the 
TCI and calculate the surface spectral function based 
on the recursive Green’s function method. It is found 
that the STCI can host mirror-protected SABSs. These 
SABSs exhibit twisted dispersion due to the surface 
Dirac cone in the normal state. In contrast to STIs, 
the position of gapless points of SABS is not fixed 
to the time reversal invariant momenta but on the 
mirror invariant lines in the surface Brillouin zone. 
Thus, the position of the gapless points depends on 
the material parameters such as chemical potential. 
As a result, the dispersion of the SABSs in the STCI 
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FIG. 1: (a)Crystal strcture of SnTe. (b) First Brillouin zone 
of a rock salt crystal and, (001) and (111) surface Brillouin 
zone. The shaded planes are (110) mirror invariant planes, 
(c) Schematic picture of the surface double Dirac cone on the 
(001) surface. 


is not a zero-energy flat one which has been seen in 
STIs, but a more complicated one. To understand the 
topological nature of this material in more detail, we 
also calculate the mirror Chern number and construct 
the phase diagram in terms of the magnitudes of the 
pairing potential and the spin-orbit interaction. We find 
that the mirror Chern numbers um takes three kinds 
of values, um = 0, —2 and —4, in STCI. Those val¬ 
ues are different from those in STIs with um = 0, — 1, — 2. 

This paper is organized as follows. In Sec. II, we re¬ 
view a tight-binding model proposed by Lent et. al. for 
IV-VI semiconductors and propose possible pair poten¬ 
tials for the tight-binding model. The topological nature 
of the TCI with possible pair potentials is discussed in 
Sec. III. In Sec. IV, we show our calculated results of 
spectral function for the (001) surface and discuss the 
evolution of the mirror-protected SABS. In Sec. V, we 
discuss the experimental situation and possibility of time- 
reversal breaking pairing. Finally, we summarize our re¬ 
sults in Sec. VI. 


II. MODEL 

Before the introduction of a model, we briefly review 
the topological crystalline insulator SnTe [5l|. In FigdJ 
we show the crystal structure (a) and, the first Brillouin 
zone and the surface Brillouin zone (b) of SnTe. The 
crystal of SnTe is the rock-salt structure, which belongs 
to the Oh point group. The band inversion occurs at four 


L points in SnTe. For the TCIs, the (001) surface is par¬ 
ticularly unique because two L points are projected to the 
same X points on the surface Brillouin zone. In this case, 
two surface Dirac cones show up at the same X points 
with different energies. The Dirac cones mix with each 
other and thus have gap in general. However, on F — X 
mirror-symmetric line, two Dirac cones remain gapless 
because each Dirac cone belongs to different mirror sub¬ 
space. A schematic picture of the Dirac cone is shown 
in Fig. [1] (c). These surface states have been observed 
by angle-resolved photoemission spectroscopy (ARPES) 
experiments [52,[53- It is noted that even if In is doped, 
SnTe still has the gapless Dirac cones [b^. 

To describe the normal state of the TCI, we use a tight- 
binding model proposed by Lent et. al. [^. This model 
Hamiltonian is a 36 x 36 matrix with nine orbitals, i.e. 5, 
two spin degrees of freedom, and two sublattices 
(anion and cation). The Hamiltonian has the following 
nearest-neighbor hopping and the on-site spin-orbit in¬ 
teraction, which is given by 
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0 
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Here, indices a and c denote anion and cation, 
respectively, and we take the basis as follows, 

(Cs,c,tTs,c,tTs,a,tTs,a,tTpx,c,tTpy,c,tTp;s,c,tTpa;,c,tTpy,c,t5 

^^ 3 z 2 _r '2 ’*^’1 T<ia;y,C,t iC'dzx,C,"[ 5C5I’^<^3^2 _,.2 jCjt? 

^dzx,a,"[^^d^2_y2 ,a,tT<i 3 ^ 2_^2 ,a,i^Cd^y ^a,i^Cdy^^a,i^Cdzx,Ci,i)' 

See the original paper for the explicit form of the 
Hamiltonian and the parameters for SnTe [6l|. The 
BdG Hamiltonian for the superconducting state is given 
by 


-HbcIG — [-HsnTe T (2) 

where represents the Pauli matrix in the Nambu 
space, ji is the chemical potential, I is an iden¬ 
tity matrix and is the pair potential. The 
basis is taken as (ca,,3,t, 

CM, Px') Py) Pz-) d^^2_j.2 ^ dy 2 ;') dzx and 

/3 = a, c. 

Up to now, almost all experiments for smerconduct¬ 
ing state in Sni_a,Ina,Te, e.g., specific heat [6214641 . //-SR 
measurements [65|, and thermal conductivity[66|, have 
suggested the fully gapped superconductivity. One of the 
candidates of the fully-gapped superconducting state is 
s-wave pair potential, which is a spin-singlet even-parity 
on-site pair potential. This pair potential is essentially 
the same as BCS state, and therefore, the energy-gap 
structure is fully gapped. Hereafter, we denote this pair 
potential as Ai. Other than the above BCS state, one 
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of the p-wave states realizes the fully gapped state. For 
p-wave pair potential in the cubic lattice, there are four 
types of time-reversal-invariant pai ring , Aiu^ Eu 
and T 2 n, as summarized in Refs. |67H7Q|. The energy gap 
of Aiu pairing is fully gapped. While the other pairings 
give point or line nodes. Thus, according to the exper¬ 
imental situation, we focus on Aiu pair potential, and 
hereafter, we denote the Aiu pair potential as A 2 . A 2 
is a spin-triplet odd-parity one, which is similar to that 
for the BW phase of ^He. In the A 2 state, the Cooper 
pair is formed between the nearest neighbor sites, i.e., 
Sn and Te sites. Following, we assume that Cooper pairs 
are formed within the p orbitals and ignore the s and 
d orbitals since low energy band is mainly composed of 
p orbitals. For simplicity, we also assume that electrons 
in Pi orbital (i = x^y^z) form Cooper pair along the i- 
direction. These simplifications do not give a qualitative 
change of the results. The explicit form of the pair po¬ 
tentials is given by 



(3) 

i=x,y,z 


\ -^ y. 

A 2^ sm—piSiax, 

(4) 

i=x,y,z 



where Si and di are the Pauli matrices of the spin and 
sublatice, respectively. In this basis, .sq [y = x^y^z)] 
denotes the spin singlet (spin triplet) and <Jo {(Jx) de¬ 
notes the intrasitre (intersite) pairing, pi is the following 
matrix: 


Pair potential 

Oh 

spin I M(iio) 

energy gap 

Ai 

A 2 

Alg 

Aiu 

singlet + + 
triplet — — 

full gap 
full gap 


TABLE I: Possible pair potentials for the STCI. Ai (A 2 ) is 
even (odd) under inversion operation. Ai (A 2 ) is even (odd) 
under the (110) mirror reflection operation. 


exploit form of Ado is given by. 


Mo 


/I 000 0 000 0\ 
0010 0 0000 
0100 0 0000 
0001 0 0000 
0000 -1 0000 
0000 0 1000 
0000 0 0100 
0000 0 0001 
Vo 0 0 0 0 0010/ 


(9) 


Since Ai has an even-parity under the (110) mirror 
operation, iLedG commutes with 


•^^10) 


•Ad(iio) 0 ^ 

0 


( 10 ) 


in the mirror-invariant plane. On the other hand, A 2 has 
an odd parity under the (110) mirror operation. There¬ 
fore, i^BdG can be block diagonalized under the diagonal 
basis of a mirror operator 


Pi 


/OO 0 00000 0\ 

0 5ix 0 0 00000 

00 000000 
0 0 O^i^ 00000 

00 0 000000 
00 0 000000 
00 0 000000 
00 0 000000 
Voo 0 000000 / 


M 


( 110 ) 




( 11 ) 


We summarize the characters of the pair potentials in 
Table H 


III. TOPOLOGICAL NATURE 


in the basis, 

i^dyz^yi^dzx^l3 1 whcrc /3 a, c and 
7 =t, The pair potentials also have the mirror reflec¬ 
tion symmetry with respect to the (110) plane. Ai (A 2 ) 
is even (odd) under the mirror reflection 


Ad(iio) AiAd|^^Q^ — Ai, 

(6) 

Ad(iio) A2Ad|^^Q^ = —A2. 

(7) 


Here, Ad (no) is the (110) mirror operator given by 


Ad 


( 110 ) 




Sy) (g) (To ( 8 ) Ado 


( 8 ) 


where Ado is a mirror operator for the orbital space. The 


In this section, we clarify the topological nature of the 
STCI, based on the recently expanded topological peri¬ 
odic table that includes crystalline symmetries in addi¬ 
tion to time-reversal and particle-hole symmetries [nl- 
[ 7 ^. First of all, we briefly review the topological nature 
of the normal state. In terms of Alt land-Zirnbauer (AZ) 
classification [iHzl, the TCI belongs to class All, since 
it has the time-reversal symmetry for fermions. However, 
the corresponding Z 2 topological number is trivial, be¬ 
cause band inversion occurs at an even number (four) of 
time-reversal invariant momenta. Thus, the time-reversal 
symmetry is not sufficient to explain the surface states in 
the TCI. Instead, the mirror reflection symmetry, which 
is specific to the crystal structure of the TCI, is responsi¬ 
ble for the stability of the surface states in the TCI. Tak¬ 
ing into account the mirror reflection symmetry, SnTe is 
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A/AsnTe 

FIG. 2: Phase diagram of superconducting topological crys¬ 
talline insulator with the odd-parity full-gapped pair potential 
A 2 as a function of the magnitude of the pair potential A and 
spin-orbit interaction A/AsnTe- Color shows the magnitude of 
the energy gap at the L point. Red dotted lines show gap 
closed line. 



classified as class All with with d = 3 and d\\ = 1 in 
Ref. [^, because the mirror operator Ad(iio) commutes 
with time-reversal operator T and obeys = — 1. 

Therefore, its topological nature is characterized by the 
mirror Chern number 


Here, n±i is the Chern number for mirror subsectors la¬ 
beled with mirror eigenvalues It has been known that 
the mirror Chern number for the TCI is —2 [5ll |. 

Next, we see the topological nature in the supercon¬ 
ducting state. In terms of the AZ classification, the 
STCI with A 2 belongs to class DHL The superconduct¬ 
ing state also possesses the mirror reflection symmetry. 
The mirror operator ot)eys (^^ 10 ))^ ~ 

commutes (anticomutes) with time reversal operator T 
(chiral operator C). Thus, this system belongs to class 
Dill with U^_ with d = 3 and dy = 1. In this case 
topological nature is characterized with two topological 
numbers, Z 0 Z. One of the integer numbers is three- 
dimensional winding number. Although the exact cal¬ 
culation of the winding number is difficult because of 
the huge Hamiltonian matrix, its parity can be easily ob¬ 
tained by using Fermi surface criteria . In the present 
case, the Fermi surfaces enclose four time reversal invari¬ 
ant momenta, i.e., L points. Thus, the winding number 
must be an even number. The other topological number 
is the mirror Chern number which can be defined on a 
two-dimensional mirror plane in the Brillouin zone. We 
calculate the mirror Chern number as a function of the 
magnitude of the spin-orbit interaction and the pair po¬ 
tential using the method proposed by Suzuki et. al [^ . 
In Fig 121 we show the topological phase diagram of the 
STCI with A 2 . The color shows the band gap at the 
L point. Here, we set fi = 0.27 eV. The mirror Chern 
number changes — 4 ^ — 2 0 with gap closing. This is 
different from the case in the STI, since the mirror Chern 
number of the STI changes —2 ^ — 1 ^ 0 j45| . 


FIG. 3: Surface spectral function for the (001) surface, (a) 
and (b) show the spectral function for the even-parity spin- 
singlet pair potential Ai along mirror symmetric F-A-F line 
and M-X-M line, respectively, (c) and (d) show the spec¬ 
tral function for the odd-parity spin-triplet pair potential A 2 
along mirror symmetric F-A-F line and M-X-M line, respec¬ 
tively. In the case of Ai, there is no Andreev bound state 
in the superconducting gap. On the other hand, there are 
mirror-protected zero-energy Andreev bound state. 


For A 2 , the BdG Hamiltonian on the mirror invariant 
plane can be divided into two eigenvectors of the mir¬ 
ror operator in Eq. m, each of which belongs to class 
D. This state is the topological crystalline supercon¬ 
ducting state defined in Ref. [s^ since each subsector 
has particle-hole symmetry. 


IV. SURFACE SPECTRAL FUNCTION 

In this section, we present our numerically calculated 
results of the surface spectral function for the (001) sur¬ 
face. To obtain the surface spectral function, we use the 



FIG. 4: Mirror separated spectral function in the case of A 2 
for (a) +z and (b) —i along the mirror symmetric line F-A-F. 
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FIG. 5: Systematic change of the surface Andreev bound states (SABSs) in the superconducting topological insulator (STI) 
(a) and the superconducting topological crystalline insulator (STCI) (b). Blue solid (red dotted) lines show the SABSs for the 
mirror +z (—z) sector. The SABSs change depending on the position of the chemical potential. Here, we assume CuxBi 2 Se 3 
and Sni-ccIuccTe for STI and STCI, respectively. For this reason, the figure of normal state for the STI (STCI) is depicted as 
the electron- (hole-) doped case. In both cases, carrier number decreases from the panel (I) to (IV). 


recursive Green’s function method using Mobius trans¬ 
formation in the semi-finite system proposed by Umerski 
[8l| . See Supplemental Material for the details of the cal¬ 
culation method. In this calculation, we set A = 0.06 eV 
in Eqs. m and The chemical potential is taken as 
/i = — 0.2 eV to fit the size of Fermi surface observed in 
ARPES measurements [g^. In this case, the conduction 
band and the TCI surface state are well separated at the 
Eermi level, and the mirror Chern number is —2. 

Eirst, we show the surface spectral function for Ai. 
Eigures. [3] (a) and (b) show the spectral function along 
the F-X-F and M-X-M lines, respectively. In the case 
of Ai, both bulk and TCI surface states have a super¬ 
conducting gap and there is no inner gap state. Next, 
in Figs. [3] (c) and (d), we show the surface spectral 
function for A 2 . On the mirror symmetric line shown in 
Fig. [3](c), there are two zero-energy SABSs correspond¬ 
ing to the mirror Chern number = —2. Different from 
the simple linear dispersion in the TCI, the dispersion of 
the surface state in the superconducting state is twisted 
since the SABSs merge into the Dirac cone in the normal 
state which exist outside the Fermi wave vector. Using 
the mirror reflection operator, we decompose the surface 
spectral function into two subsectors with different eigen¬ 
value ±i as shown in Fig. 01 As we mentioned in Sec. 
M the Chern number of each subsector is =f 2. There¬ 
fore, there are two chiral SABSs in each subsector. There 
are also zero-energy SABSs along M-X-M line [see Fig. 
[3Kd)]. These zero energy SABSs can be interpreted by 
zero-dimensional topological number defined by combin¬ 
ing with Ad( 110 ) and C [^. These zero-energy SABSs 
disappear if we enlarge the magnitude of A like A = 0.1 
eV. 

We now turn to the detail analysis on the mirror pro¬ 


tected SABSs which appears in the case of A 2 . First, we 
briefly review the systematic change of the time-reversal 
symmetry (TRS) protected SABSs in the STI such as 
Cua^Bi 2 Se 3 [43-[46|. The dispersion of the SABSs in the 
STI changes mainly with the chemical potential. If the 
Fermi level is much higher than the energy where surface 
state is merged into bulk band, the shape of the SABSs is 
normal cone shape similar to that of the ^He BW phase, 
as can be seen in Fig. [5] (a) (I). On the other hand, as 
lower the Fermi level, the SABSs start twisting as shown 
in Fig. [5] (a) (III) and, finally, the SABSs become like 
Fig. [ 5 ] (a) (IV). This is regarded as the Lifshitz transition 
of the surface state. At the transition point of this Lif¬ 
shitz transition, the dispersion becomes flat like as shown 
in Fig. [5] (a) (II) since the position of zero energy SABS 
is pinned at time reversal invariant momenta, i. e. k = 0. 



FIG. 6: Surface density of state in the case of A 2 for the (001) 
surface. There are four peaks at E/A = ±0.49 and ±0.16 















































6 




FIG. 8: Surface spectral function in the case of A 2 + iAi 
pairing. We introduce Ai and A 2 at the same ratio. 


V. INTERPRETATION OF EXPERIMENTS 
AND DISCUSSION 


FIG. 7: Surface spectral function for the (111) surface along 
mirror-symmetric line F-M. (a) and (b) show the sur¬ 
face spectral function of the normal state for Sn- and Te- 
terminated surface, respectively, (c) and (d) show the spectral 
function of superconducting state for Sn- and Te-terminated 
surfaces, respectively. Here, we set the chemical potential 
/i = —0.2 eV, which is the same as the calculation for the 
(001) surface. 


This flat-like band can explain the ZB CP in the point- 
contact experiments even if the bulk energy spectrum is 
fully gapped 

Next, we reveal the structural transition of the SABSs 
in the STCI. If the Fermi level is much lower than the 
energy where the TCI surface state is merged into the 
bulk band, there are two Dirac cones on the mirror 
plane as shown in Fig. [5] (b) (I). On the other hand, in 
the same manner as the STI, if the chemical potential 
is located on the energy where surface and bulk states 
are well separated, the dispersion of the SABSs in the 
STCI is twisted as shown in Fig. [5fb)(IV). However, 
different from the SABSs in the STI, this change of 
the dispersion does not involve the Lifshitz transition. 
Thus, the surface state of the STCI does not have to 
host zero-energy flat dispersion between the twisting (I) 
and non-twisting (IV) cases as shown in Fig. [51(b) (II) 
and (III). This can be understood by the difference of 
the symmetry which protects the zero energy SABSs. 
Different from the case of the STI, the zero-energy 
SABSs in the STCI are not protected by the TRS but 
the mirror symmetry. Therefore, the position of the 
zero-energy SABSs in the surface Brillouin zone can 
move along the mirror-symmetric line with doping. For 
this reason, the mirror-protected SABSs in the STCI can 
twist without presenting the flat-like band at zero energy. 


We discuss the experimental situation for the tunneling 
spectroscopy. First, in Fig. |6l we show our numerically 
calculated results of normalized surface density of state 
ps{E)/pn{^) for the (001) surface in the case of A 2 . Here, 
we set A = 0.02 eV. As can be seen from FiglHl there 
are four peaks at E'/A = ±0.49 and ±0.16. Therefore, 
it is reasonable that four peaks appear in conductance 
measurement with low transmissivity realized in scan¬ 
ning tunneling spectroscopy (STS). On the other hand, 
in the point-contact measurements, it is likely to include 
the components of the surface with different orientations 
since the surface of crystal may contain many steps in the 
contact region. Thus, we calculate the surface spectral 
function for the (111) surface as an example of different 
surfaces. In the case of the (111) surface, L points are 
projected to either F point or M points as shown in Fig. 
mb). Moreover, there are two distinct situations, i.e., Sn- 
terminated surface and Te-terminated surface. In Fig. 0 
(a) [(b)], the spectral function for Sn- (Te-) terminated 
surface in the normal state is presented. In the case of 
the Sn-terminated surface, Dirac points are located on 
the high symmetric points F and M. On the other hand, 
in the case of Te-terminated surface, Dirac points appear 
on the mirror symmetric line F-M. In FigH (c)((d)), we 
show the Sn(Te)-terminated surface spectral function in 
the superconducting state for A 2 . As can be seen from 
the figures, the SABSs for the (111) surface is completely 
different from those for the (001) surface. In the case of 
Sn-terminated surface, two cone-shape surface states ap¬ 
pear at the F and M points on the surface Brillouin zone 
corresponding to mirror Chern number um = —2. In this 
case, the change of the dispersion of the SABS is almost 
the same as that of the STI. In the case of Te-terminated 
surface, as you can see in Figl7|(d), the surface state at M 
point is flat like and there are zero-energy surface states 
between the F and M points since the valley surface state 
remains due to the existence of the mirror reflection sym¬ 
metry. 

From the above calculation, it is reasonable to conclude 
that the bias voltage dependence of conductance for the 
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(001) surface can be changed from the four peaks spectra 
if the component of the (111) surface is included. In ad¬ 
dition, in the normal state, it has been revealed that the 
edge state of (001) film dramatically changes depending 
on the number of layers [HI . It can be expected that the 
similar situation may occur in the superconducting state 
due to the existence of the steps on the surface. There¬ 
fore, further calculations are necessary to interpret the 
obtained point-contact experiments [HI- 

Before closing this section, we briefly refer to the case 
of A 2 + iAi, i.e., p-\-is pairing which breaks time-reversal 
symmetry. Novak et. al. have indicated the realization 
of p -1- impairing in Sni-a^In^^Te by point-contact exper¬ 
iments [HI- Theoretically, Goswami and Roy have sug¬ 
gested the possibility of the p + is pairing in doped topo¬ 
logical materials with an effective model [HI . Motivated 
by the studies, we calculate the spectral function for the 
A 2 + iAi state in the lattice model. The calculated re¬ 
sults are shown in Fig. [51 Consequently, we find that 
even in the lattice model, the system has gapped SABSs 
since A 2 +iAi pairing breaks the (110) mirror symmetry. 
Our model is useful for further discussion of anomalous 
thermal Hall effect in Sni-^^In^^Te. 


VI. CONCLUSION 

In this paper, we have theoretically studied the SnTe 
class TCI in the superconducting state. We have intro¬ 
duced the possible fully-gapped pair potentials to the 
TCI and calculated the surface spectral function using 
the recursive Green’s function method. We have found 
that the STCI hosts the mirror-protected SABSs when 
odd-parity pair potential A 2 is realized. The mirror- 


protected SABSs for the (001) surface appear only when 
we use the model that describes the entire Brillouin zone. 
This result suggests the importance of considering the 
overlap of the Fermi surface to reveal the SABSs in gen¬ 
eral. We have also shown the systematic change of the 
dispersion of the mirror-protected SABSs. It has been re¬ 
vealed that the dispersion of the mirror-protected SABSs 
is not necessary to become flat at zero energy when it 
twists, which is different from the SABSs in the STL 
Moreover, we find that four peaks appear in the SDOS, 
which can be detected by STS measurements. 

In analogy with the mirror-protected Dirac cone in the 
normal state, the mirror-protected SABSs can be tuned 
with the electronic field or strain. This means that some 
physical properties, e.g., the Josephson current can be 
controlled by these effects. We expect that further study 
of the mirror-protected SABS will lead to the develop¬ 
ment of superconducting electronics. 
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